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Double-pair emission from typc-II parametric down conversion results in a highly entangled 4- 
photon state. Due to interference, which is similar to bunching from thermal emission, this state is 
not simply a product of two pairs. The observation of this state can be achieved by splitting the 
two emission modes at beam splitters and subsequent detection of a photon in each output. Here 
we describe the features of this state and give a Bell theorem for a 4-photon test of local realistic 
hidden variable theories. 

PACS Numbers: 3.65 Bz, 3.67 -a, 42.50 Ar 

Parametric down conversion has proven to be the best source of entangled photon pairs so far in an ever increasing 
^ ■ number of experiments on the foundations of quantum mechanics ||l[ and in the new field of quantum communication. 
Experimental realizations of concepts like entanglement based quantum cryptography j2| , quantum teleportation j^] 
and its variations |^ demonstrated the usability of this source. New proposals for quantum communication schemes 
' 1^ and, of course, for improved tests of local hidden variable theories initiated the quest for entangled multi-photon 
states. Interference of photons generated by independent down conversion processes enabled the first demonstration 
of a three-photon Greenberger-Horne-Zeilinger (GHZ)-argument j6j and, quite recently, even the observation of a 
> ; four-photon GHZ-state §. 

In this report we show that four photon entanglement can be obtained directly from type-II parametric down 
\ conversion. Instead of sophisticated but fragile interferometric set-ups, we utilize bosonic interference in a double-pair 
' emission process. This effect causes strong correlations between measurement results of the 4 photons and renders 
■ type-II down conversion a valuable tool for new multi party quantum communication schemes. The analysis of the 
^-H entanglement inherent in the four photon emission leads us to a new form of inequality distinguishing local hidden 
variable theories from quantum mechanics, and demonstrates its potentiality for experiments on the foundations of 
quantum mechanics. 

In type-II parametric down conversion Q multiple emission events during a single pump pulse lead to the following 
I ' state 

^' . . . 

where Z is a normalization constant, a is proportional to the pulse amplitude, and where ay is the creation operator 
of a photon with vertical polarization in mode a, etc. (Fig. 1). We are interested only in 4 photon effects, i.e. the 
emission of two pairs. Then only the term in (m) proportional to 
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{a*Hb*y + a*yb*Hf\0), (2) 
is relevant. The particle interpretation of this term can be obtained by its expansion 

{a*j!b*y^ + a*^b*^ + 2a*ya*Hb*yb*H) \ 0) , (3) 
and is given by the following superposition of photon number states 

I 2Ha, 2Vfc) + I 2Va, 2Hb) + \ IHa, IK, IHb, IVb) , (4) 

where e.g. 2Ha means 2 H polarized photons in the beam a. 

One should stress here that this type of description is valid only for down conversion emissions, which is detected 
behind filters endowed with a frequency band, which is narrower than that of the pumping fields ||. If a wide band 
down-conversion is accepted then such a state is effective only if counts at the detectors are treated as coincidences, 
when they occur within time windows narrower than the inverse of the bandwidth of the radiation If such 

conditions are not met, then the four photon events are essentially emissions of two independent, entangled pairs, 
with the entanglement existing only within each pair. 

Let us pass the four photon state via two polarization independent 50 — 50 beam splitters. For simplicity we assume 
that at the beam splitters a is transformed into "^(a + a') and b into -^[b + b'), with prime denoting the reflected 
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beam. One can expand the expression (Q), and then extract only those terms that lead to 4 photon coincidence 
behind the two beam splitters, i.e. only those terms for which there is one photon in each of the beams. The resulting 
component of the full state is given by 

[4 {a*Ha';jb*yb'^ + a*ya'^b*Hh'}i) + 2(a^,a^ + a*ya'^){b*yb'^ + b*^^^)] \ 0) . (5) 

The first term represents a 4-photon GHZ state, whereas the second one is a product state of two EPR-Bohm states 
(the I 4'"'") Bell states in polarizations H and V). This, after the normalization, can be symbolically written as 

V273 I GHZ)^^,,,, + 7173 I EPR)^^, I EPR),,, . (6) 



'bb' y ^<.laa' I ^'■/bb' ■ 

For additional simplicity of the presentation we also rotate the polarizations in the beams a and a' by 90° . Thus now 
our initial state is given by with the GHZ state in its standard form, resulting in the state 

yi73 ( I VVVV)^^,,,, + I HHHH)^^,,,, 
+ i (I HVHV)^^,,,, + I HVVH)^^.,,, + I VHHV)^^,,,. + \ V HV H) ^^,,,,)) . (7) 

In order to demonstrate the entanglement of this state let us analyze polarization correlation measurements involving 
all four exit ports of the beam splitters, where the actual observables to be measured are elliptic polarizations with 
main axis of the polarization ellipse at 45°. Such observables are of dichotomic nature, i.e. endowed with two valued 
spectrum k = +1,-1, and are defined for each spatial propagation mode x = a, a', 5, 6' by their eigenstates 

s/Tj2\V),^ + ke~'*^^\H)^ = \kA.)- (8) 

The probability amplitudes for the results k,l,m,n = ±1 at the detector stations in the beams a,a',b,b', under 
local phase settings (pa, (pa' , 4>b, 4>b' , respectively, are given by 
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4a/3 



1 + klmn E ^ + i (fc e"^" + / e*»' ) (m e"^" + n e*"^''' ) 



(9) 



where 4> stands for the sum of all local phases. Therefore the probability to get a particular set of results {k,l,m,n) 
is given by 

P{k, I, m, n\(f>a, (pa' , 06, (pb') 
= + klmncos^(j}) 

+ + klcos {(pa " 0a'))(l + mncos{(j)b - (pb')) 
+ lRe(^{l + klmne'T.'t'){ke"t'- +le"^-'){me"^'' +ne'l'»')^ . (10) 

The last term is written in the form of a real part of a complex function to shorten the expression. 
The correlation function is defined as the mean value of the product of the four local results 

E{(pa,(pa',(l>b,(pb') = Efc=±l Ei=±l Em=±l E„=±l " (H) 

Its explicit form for the considered process is given by 

E{(pa,(pa',(pb,(pb') = I COSX;^ + I COS (0a - (pa') COS {(pb " Pb') ■ (12) 

Only the first two terms of the probabilities ([lo| ) contribute to the correlation function, and the function is itself a 
weighted sum of the GHZ correlation function (the first term) and a product of two EPR-Bell correlation functions. 
The last term in (^o|) gives a zero contribution to the correlation function, because sums like EfcZmn vanish. 

The correlation function (^2|) for the process has a more complicated form than in the usual cases for GHZ-type 
states, but the strong correlations for numerous phase settings clearly indicate incompatibility with local realistic 
theories. When inserted into Mermin-type Bell inequalities for four particle systems, the violation is not too 
impressive, even for optimal sets of local phases. However, here we present a reasoning, involving Bell inequalities of 
a new type , giving stronger inequalities for distinguishing the validity of the different theories in a four photon 
experiment. 

In a local hidden variable theory a correlation function has to be modeled by a construction of the following form 
(see e.g. M) 
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ELHv{(l>a,4>a',4'b,4'b') = J d\p{X)Ia{4>a, ^^la' {4>a' A)h{<t}bA)h' [4>b' , , (13) 

where A represents an arbitrary set of values of local hidden variables, p{X) their probabilistic distribution, and 
Ixiipx, A) = ±1 (x = a, a' , 6, b') represent the predetermined values of the measurements. Their values depend on the 
set of hidden variables and on the value of the local phase settings. 

We start with allowing each observer of beam x (= a, a', b, b') to choose, just like in the standard cases of the Bell 
and GHZ theorems ||l^, between two values (pl. and i/)^ of the local phase settings. 

The formula for the LHV correlation function for the chosen settings is given by 

ELHv{ra,<l>l',cl,l, <j>l,)^J dXp{X)Ia{ra,^)Ia' {cl>l, , X) h {cl>l , X)!,, [cj^l, , X) (14) 

with p, q,r,s = 1, 2. It is important to stress that one must consider arbitrary LHV correlation functions. The only 
constraint being their structure given by (p^). 

One can treat the full set of the LHV predictions as a four index tensor ElhVi with the indices p,q,r,s — 1,2, 
built out of the tensorial products of two dimensional real vectors = ( Jj, (</>;[,, A) , /a; {(p1 , A) ) , which represent the two 
possible results of a given observer for the given value of the hidden variable: 

Elhv = J dXp{X)v^ (g) v^, (g)v^(g) v^, (15) 

The actual values of the components of the two dimensional vectors {Ix{(f>l:, A), Ixifft^ •^)) can be equal to only either 
(1,1), or (1,-1), or (—1,-1), or finally (—1,1). Let us denote these four possible vectors by v^^. with j — 1,2,3,4, 
respectively. Thus, the LHV correlation function (tensor) can be simplified to a discrete sum over hidden probabilities 
Pk,i,m,n of the tensorial products of all possible measurement results. 

Elhv = Ek.Lm.n=K..A Pk,l,m,n ® v^, ® ® vg, (16) 

A further simplification of the tensor is possible since (—1, —1) = —(1, 1) and (—1, 1) = —(1, —1), or in other words 
v^+^ — — v^. The tensorial products (g) v^, (g) v™ (g) v{J, with k, l,m,n = 1,2 form a complete orthogonal (product) 
basis in the (real) Hilbert space of tensors E? ® <g) ® E? . One can thus rewrite the expansion (|l^ ) so that it 
becomes an expansion in terms of the aforementioned basis 

Elhv = Efc,i,™,„=i,2 Cfc,i,m,„ (g v^, v™ ® v^, . (17) 

The relation between the coefficients in and the probabihties of (^6|) is given by 

CA:,/,m,n Vk,l.rn;n Pk-\-2,l.7n,n Pk.l-\-2,m.7i ■■■ 

+ ••• - Pk+2,l+2,m+2,n + ■■■ + Pk+2,l+2jn+2,m+2 ■ (18) 

The expansion coefficients are of course unique, and since X^fc z m n=i 4Pk,i,m,n = 1 they satisfy the following 
inequality 

Efc,i,m,„=l,2 |c/c,i,m,n| < 1 • (19) 

This inequality is a necessary condition for the local realistic description to hold, and thus gives the handle for 
evaluating the validity of this class of theories. It should be stressed that one can also show that the inequality ( |l9| ) 
is also a sufficient condition (see [Q). 

To compare the structure of the possible LHV correlation functions with our quantum one (p^), let us, in order to 
simplify the analysis, choose specific values for the local phase settings. First, the observer of beam a will be allowed 
the choice between cj)]^ — and i/)^ = tt/2. The other observers {y = a',b,b') can choose between (py = — 7r/4 and 
(j)y = 7r/4. Next, one can expand the quantum function ( p^ ) into a sum of products of sine and cosine functions of 
single phases 

E{<j)a,4'a',4>b, <Pb') 
= CaCa'CbCb' + SaSa'SbSb' - \{SaSa'ChCv + CaCa' SbSy) 

- 'i{SaCa'SbCb' + CaSa'CbSb' + CaSa'SbCb' + SaCa'CbSb') , (20) 



where Sx — smcpx and Cx = cos . For each fixed set of four local settings one can calculate the specific value of the 



quantum correlation fimction i?((/)^, 0^, , (/)^, 0^ ). We notice, that for the specific phase settings given above one has 
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and 



whereas 



and 



(sin(<^i),sin(</)2)) = -i=(l,-l), 



(cos (01), cos - (1, 0) = 1) + -1) 



(sin(0i),sin(<^2)) ^ (0^1) ^ hi^^) _ 1(1, 



Therefore the quantum predictions can be arranged to form a tensor, too, and it is easy to write down its expansion 
in the product basis (the same basis as in eq. (|l7|)) 

E = Efc,i.™,„=l,2 9fe,i,™,nvk ^ ^1^, ^ ^ ^ (21) 



The actual values of the expansion coefficients qk,i,m,n can be straightforwardly obtained from (20). However, note 
that for the specific set of angles chosen above one has 

Sfc,i,m,n=l,2 \lk,l,m,n\ — > 1 ■ (22) 

Keeping in mind that the expansion in terms of basis vectors is unique, the quantum correlation function out of which 
the tensor E is built thus violates the necessarycondition for local realism (p^). 

The quantum correlation function satisfies (^^ only if it is multiplied by a scaling factor v equal or smaller than 

~ 53%, in other words, if one replaces it by 

E'i(t)a, 4>a',<t)b, 4>b') = vE{(j)a,4'a',4>b, 06') ■ (23) 

In an interferometry experiment this scaling parameter is directly related to the visibility (contrast) of the interference 
pattern. Visibilities lower than one can be interpreted as arising due to some noise contribution to the state. If one 
considers mixed states for the system of the type 

Pv = {'^ - v)pnoise + {ipl , (24) 

where Pnoise = Jq^ represents completely uncorrelated noise contribution, and | ip) stands for our pure state (^), then 
the aforementioned critical v gives the threshold beyond which no LHV model can resemble the quantum predictions. 

In conclusion, parametric type-II down conversion not only produces entangled photon pairs, but also highly 
entangled four photon states. The observation of these states is experimentally much easier to achieve than for GHZ- 
type states. Here, the full set of probabilities of possible LHV-predictions is compared with the quantum predictions 
resulting in a significant distinction of the theories. 

The other interesting feature of the considered process is that for a number of specific settings one obtains perfect 
four photon correlations. E.g., for all local phases equal to zero the correlation function is equal to 1, whereas for 
4'b = 0b', 0a — 0a' = and 0a' + 06 = it is cqual to —1. This directly enables one to transfer the standard 
protocols for entanglement based quantum cryptography ||l5|] to the four photon case making multi party quantum 
key distribution and quantum communication complexity schemes feasible. 

MZ was supported by University of Gdansk grant BW-5400-5-0032-0. and the Erwin Schrodinger International 
Institute for Mathematical Physics, Austria. This work was supported by the EU-Project QuComm (IST-1999- 
10033). 
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FIG. 1. Experimental set up to demonstrate the entanglement inherent in 4-photon emission from type-II paramteric down 
conversion. (A/2-half wave plate to flip polarization, BS-beam splitter, and P^^ represents polarization analysis corresponding 
to the phase angle (pa, etc.) 
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